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Maleev-like boson transformation. The ground-state properties and the triplet excitation spectra are 
obtained as functions of the coupling parameter. Secondly, based on the pseudo-spin operators, a 
microscopic treatment is formulated within the framework of the powerful, systematic coupled-cluster 
method (CCM). Comparison between various approximations is made. The advantage of the CCM for 
the purposes of systematic improvement is emphasized. 
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I. Introduction 



Spontaneous dimerization of theoretical spin-lattice models was perhaps first discovered by Ma- 
jumdar and Ghosh [f] in 1969. They found that for the one-dimensional (ID) spin-i system with 
nearest-neighbour and next-nearest-neighbour couplings, the perfect dimer state, in which every two 
adjacent atoms form a spin singlet valence-bond, is the exact ground state of the Hamiltonian at a 
particular coupling. Obviously, the translational symmetry is spontaneously broken in the dimerized 
system, and the corresponding ground state of the ID model is doubly degenerate. 

A spin-1 chain or other higher-order spin chains can also exhibit dimerized valence-bond structures 
in their ground states [2]. Very recently, a series of ID SU{n) spin-s (n — 2s + 1) antifcrromagnetic 
models have been solved by Bethe ansatz [3] , and the corresponding dimerization order parameters have 
been exactly calculated by the author [4]. In addition to the dimerized spin chains, it is possible that 
some ID spin systems (with integer spin quantum numbers) favour the trimer configuration which is 
produced by a sequence of spin-singlet states formed from every three adjacent spins. Furthermore, 
dimerization or trimcrization may also occur in two or higher-order spatial dimensions. Clearly, a 
dimerized or trimerized spin system can be viewed as a kind of solid in which the corresponding simple 
valence-bonds are localized and the translational symmetry is broken. The perfect dimer or trimer state 
is not in general the exact ground state of a given Hamiltonian, but for some systems the ground state 
may still possess a nonzero solid-like dimer or trimer long-range order and hence show the characteristics 
of a quantum solid. (It therefore seems more appropriate to call these solid-like systems collectively as 
'valence-bond lattices' [5].) Due to quantum correlations, one expects that the corresponding long-range 
order of those quantum solids will be reduced or vanish at certain coupling strengths. A good example 
is provided by the ground state of the spin-1 SU(3) antiferromagnetic chain [4], where the dimerization 
order parameter is reduced to 42%. 

One is quite familiar with phonons in ordinary atomic solids and magnons in ferromagnets or an- 
tiferromagnets. In particular, the spin-wave theory of Anderson [6] provides a simple and excellent 
description of the spin correlations in the ground and low-lying excited states for a number of antiferro- 
magnets. Similar approximations have also been developed for some dimerized systems. Parkinson [7] 
formulated a spin- wave approximation based on a spin-i dimer state. He focused on the spin-i Heisen- 
berg model and employed the method of equation-of-motions. A triplet excitation spectrum of v^sinfc 
was obtained. This compares well with the exact result of triplet spectrum of -| sin k. Chubukov [8] later 
provided a similar spin-wave theory by using a transformation similar to that of Holstein-Primakoff to 
study specifically the dimerization of the spin-1 chains with Hciscnberg and biquadratic exchanges. Read 
and Sachdev [9] investigated the dimerization problem within the framework of Schwinger boson field 
theory. Very recently, the author [10] has extended Parkinson's theory to discuss possible trimerization 
of a ID isotropic spin-1 system. The concept of dimerization has also been extended to the 2D spin 
systems, e.g., the J1-J2 Hcisenberg model on the square lattice [11], etc. 

In this article, I intend to investigate the dimerization of spin systems by a systematic, microscopic 
approach. Because of its simplicity, I focus on the ID spin-i Heisenberg chains with nearest-neighbour 
and next-ncarcst-neighbour couplings, for which the dimer state is the ground state at a particular 
coupling [1]. Following Parkinson [7], I study the dimerization in terms of a proper set of composite 
operators (pseudo-spin operators). I extend and reformulate Parkinson's theory so that the ground-state 
properties as well as the excitations can be investigated. Firstly, by employing a Dyson-Maleev-like boson 
transformation [12] for the pseudo-spin operators, a spin-wave theory is developed. The advantage of 
using Dyson-Maleev transformation lies in the fact that Hamiltonian can be expressed in a compact 
form as a finite order polynomial of boson operators, rather than as an infinite series which is the case 
when the Holstein-Primakoff-like transformation is used (e.g. Ref. [6]). Secondly, I apply the powerful, 
systematic, microscopic many-body theory of the coupled-cluster method (CCM) [13] based on the 
pseudo-spin operators of the frustrated spin-i chain. A systematic approximation scheme within the 
CCM approach is developed for the ground state. The CCM has recently been successfully applied to 
the various spin systems with Ising-like long-range order [14] or with planar long-range order [15]. The 
excellent results produced by the CCM approximations, particularly for the spin systems with Ising-like 
long-range order, is the main motivation for the current CCM approach to the dimerization problems. 

It should be emphasized that the dimerized or trimerized states are not merely mathematical ar- 
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tifacts. In fact, the ID frustrated spin-^ system has recently been shown to be relevant to the recon- 
struction of fee metal surfaces at finite temperature, with the spin dimcrized phase corresponding to the 
disordered flat surfaces [16]. 

The plan of this paper is as follows. In Sec. II the composite operators in a matrix representation are 
introduced and the corresponding boson transformations are given. Sec. Ill is devoted to the spin-wave 
approximation for the ID spin-i model. The ground-state energy and excitation spectra are obtained as 
functions of the coupling constants. In Sec. IV I describe in detail the microscopic CCM for the ground 
state of the dimerized spin system. The results from the CCM approximation are compared with the 
spin-wave theory. I conclude this article by a general discussion in Sec. V. 



II. Pseudo-Spin Operators and Their Bosonizations 



I first consider a two-spin system, each with spin i. For completeness, some of the analysis by 
Parkinson [7] is repeated here. Clearly, there are four states for such a two-spin system. If notations |f) 
and | J,) are used to represent spin up and down states respectively, the singlet state can be written as 

|o) = -±=(|n)-|IT», (2.1) 

and the triplet states with sf otal (= sf + s|) = 1,0,-1 are given by respectively 

|i) = ITT), |2> = ^(|U) + UT», 13} = |U). (2.2) 

Following Parkinson [7], I employ a matrix representation. Each of the four states of Eqs. (2.1) and 
(2.2) is then represented by a column matrix with a single nonzero element. One can then introduce 
operators Aij as having only a single non-zero element of a (4 x 4) matrix, namely (i'\Aij\j') — 5w6jj>. 
For example, A\q (A 3 o) is an operator which increases (decreases) s^ otal by one unit, while A20 leaves 
it unchanged. Their hermitian conjugates (i.e., transpose matrices) have the opposite effects. Together 
with other operators, these sixteen operators form a complete set for the spin pair and any operator of 
the pair can be written as a linear combination of these sixteen. 

For a pair of spins, each of which has spin greater than one half, similar operators can be defined. 
For example, there are nine states for a two-spin system each with spin 1, and hence there are eighty-one 
(9 x 9) A^ operators which form a complete set for a pair of spin-1 atoms. For three-spin system each 
with spin 1, the dimensionality of the matrix is twenty-seven [10]. 

It is worth pointing out that these A^ operators are nonlinear in terms of the original single spin 
operators, for example, 

A 00 = |-s 1 -s 2 , A 01 = —j=(4 - 4)(4 + 4), (2-3) 

etc. It is in this sense that I have referred to these A^ as composite operators [10]. (Notice that Aoq is 
the usual spin-singlet projection operator.) Furthermore, it is easy to prove that they obey the following 
pseudo-spin algebra, 

[Aij, A kt ] = A u Sj k - A k j8 u . (2.4) 

Therefore, I also refer to them as pseudo-spin operators. My assumption in this paper is that it is more 
natural to study the dimerization in terms of these composite operators rather than the original single 
spin operators. The spin-wave theory in Sec. Ill and the CCM approximations in Sec. IV are developed 
based on this assumption. 

In a straightforward manner, one can express the single spin operators in terms of A^ operators [7] . 
They are given by [17] 

s{ = \( A 02 + A 20 + An - A 33 ), s z 2 = ^(-A) 2 - A 20 + An - A 33 ), (2.5a) 

st = -^=(A 30 - A 01 +A 21 +A 32 ), s% = -^={A i- A 30 + A 21 +A 32 ), (2.56) 

Si = -^=(A 03 - A 10 + A 12 +A 23 ), S2 - = ^=(^10 - A m + A 12 + A 23 ), (2.5c) 
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Recognizing that Aij obeys the pseudospin algebra of Eq. (2.4), one can make the following Dyson-Maleev 
transformation [12], 



Am 


= 1 — afai 


- d\d 2 — a 3 ; 






= a-pAoo, 


A@p — dp] 


(2.6) 


A 


a p dp, 


Apq dp dq, 





where p,q= 1, 2, 3, and d p , a+ are three sets of boson operators, obeying the usual boson commutation 

[dp, d+] = 1, (2.7) 

and with all other commutators yielding zero. 

By definition the singlet state |0) of Eq. (2.1) is the vacuum state of the bosons, namely, 

Op|0) = 0, p= 1,2,3. (2.8) 

The physical states correspond to the vacuum state |0) and the three states with only one boson excited. 
Furthermore, as the matrix elements between physical and unphysical subspaces are equal to zero, the 
transformation given by Eqs. (2.6) is exact at zero temperature just as in the case of the conventional 
spin- wave theory [6]. 

A general spin-i Hamiltonian can be expressed in terms of A^ operators by Eqs. (2.5) and then in 
terms of those three sets of boson operators according to Eqs. (2.6). In Sec. Ill I consider a spin- wave 
theory for the frustrated ID model using these pseudo-spin operators and their bosonizations. In Sec. IV 
I develop a microscopic formalism within the framework of the CCM also based on these pseudo-spin 
operators. 

III. Spin- Wave Theory 

The ID spin-i isotropic model with nearest-neighbour and next-nearest-neighbour couplings is 
described by the Hamiltonian, 

N 

H = ^2(si ■ S;+i + Jsi ■ s ;+2 ), (3.1) 

1=1 

where J is the coupling constant and TV is the total number of spins. I use the periodic boundary 
condition and choose even N for convenience. I have also taken the lattice spacing to be unity. At 
J = \, the Hamiltonian of Eq. (3.1) becomes the well-known Majumdar-Ghosh model [1], which has 
two degenerate dimer ground-states, with one of them given by 

N/2 

\D) = I} |0>2r-l,2r. ( 3 - 2 ) 
r=l 

where the notation |0)^ . represents the singlet-paired state of Eq. (2.1). This dimer state |D) is shown 
graphically in Fig. 1. After choosing the dimer indices as shown in Fig. 1, the Hamiltonian can be 
written as 

N/2 

H = J^[si(r) • s 2 (r) + s 2 (r) ■ si(r + 1) + Jsi(r) ■ si(r + 1) + Js 2 (r) • s 2 (r + 1)] . (3.3) 

r=l 

One can then express H in terms of the composite operators A^ by Eqs. (2.5) as 

N/2 

H = Y J {H r + - A H r , r+1 + -JH' rtr+l ), (3.4) 



with 



1 

ti r = 



Hr=\~ A r 00 , (3.5a) 



(3.56) 
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t , . = / \r I ^r- _ at _ at \i Ar+1 I Ar+1 I 4r+l _ Ar+1\ 
'r\r+l — I ^02 + ^11 ^20 /i 33A /i 02 + -^11 +^20 ^33 ) 

+ (^01 + ^21 — ^30 + y ^32)(^03 fl — ^lO 1 + 1 + -^23 *) 

+ ( — ^03 + A-iQ + + -AjaX - Af>l 1 + ^l" 1 + ^30 1 + ^S^ 1 )' 

^r,r+l = (^02 + ^oX^-oJ + ^o" 1 ) + (^11 — -^33X"^H~ — ^33 

+ (A r m Al )(A[+i - A^ 1 ) + (A r 21 + ^ 2 ) K+ 1 + A^ 1 ) 

+ (^io — ^03) (^01" — ^30 ) + (^12 + ^23) {A r 2 i + )• (3.5c) 

By Eqs. (2.6), one can further express H in terms of the three sets of boson operators. For clarity, 
I use different notations for these three sets of bosons, 

a = a\, a + = a+; b = a 3 , b + = ajj~; c=a 2 , c + = a\ . (3.6) 

Now the Hamiltonian of Eq. (3.4) can be written as 

H = H Q + V, (3.7) 

where Hq contains only the quadratic terms and a constant, 

3 N/2 r i 

H = —=N + a+a r + b+b r + c+c r + -(2 J - 1) [(a r - 6+)(a+ +1 - & r+1 ) 

8 -=i l (3.8) 

+ (a+ - 6 r )(a r+ i - 6+ +1 ) + (c r + c+)(c r+ i + c+ +1 )] |, 

and V contains higher-order terms up to the sixth, 

V = V 3 + Vi + V 5 + V 6 . (3.9) 

There is a close analogy between the bosonizations in the present case and in Anderson's spin- wave 
theory [6], where the Neel state is the vacuum state for the two sets of boson operators. Clearly, A o 
corresponds to s z in the conventional spin-wave theory, while A n0 (A 0n ) corresponds to s + (s~), etc. 
Notice that there are three sets of independent boson operators in the present case, but there are only 
two in the conventional spin-wave theory. This is because the symmetry-broken vacuum state (the 
Neel state) in the conventional spin- wave theory is in the subspace of zero Sj otal (= J2i s i) while the 
symmetry-broken vacuum state (the dimer valence-bond state) in the present case is in the subspace of 
zero vector s to tai (= J2i s 0- ^ i s already clear that one should expect a triplet excitation of spin 1 for 
the present dimerized system. 

After Fourier transformations for the three sets of boson operators, 



rjN/2 rjN/2 

flfe = ViV^ exp(-2z/cr)a r , b k = J — ^ exp(-2ifcr)& r , 



r=l 
-N/1 



Ck 



r—1 



(3.10) 



(the factor of 2 in the exponential functions is due to the double spacings in the dimer index r,) H of 
Eqs. (3.8) can be diagonalized by the usual Bogoliubov transformations, 

a k — cosh9 k a k - sinh6» fe /3+ fe , b k = - sinh6 k at k + coshO k (3 k , 

. (3.11) 
c fe = cosh^ fe 7 fc + sinh6» fe 7_ fe , 
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where Ok is given by 



(3.12) 



The diagonalizcd Hamiltonian H can be simply written as 



H = ^2u k (ala k + l3+/3 k + 7+ 7fe ) + E , 



(3.13) 



k 



where the triplet spectrum is given by 



Vl - (1 - 2 J) cos2fc 



(3.14) 



which agrees with Parkinson [7] at J = 0, 



and where Eq is defined as 



Eq 
~W 



|^[Vl-(l-2J)cos2fc-l]-| 



(3.15) 



In Eqs. (3.13) and (3.15) the summation over k is defined as 




dk. 



(3.16) 



It is well known [18] that the ground state |<j>o) of the quadratic Hamiltonian H is given by the 
two-body form as 



where 9 k is determined by Eq. (3.12) and \D) is the boson vacuum state of Eq. (3.2). 

The ground-state energy E g within this spin-wave theory is given by the expectation value of the 
full Hamiltonian of Eq. (3.7) with respect to |$o)- Clearly, the odd-body terms yield zero, and one has 
Eg = Eq + (V4) + (Via), where (V4) and (Vq) can be calculated by Wick's theorem. If one ignores V4 
and Vq which represent spin-wave interactions, the ground-state energy is then approximated by Eq of 
Eq. (3.15). This is the result shown in Fig. 2, where I also include for comparison the numerical results 
of Tonegawa and Harada [19], obtained by extrapolating the finite-size calculations for J < 1/2, and 
the exact results by Parkinson [20] of the N — 20 system for J > 1/2. Notice that at J = 1/2 (the 
Majumdar- Ghosh point), Eq. (3.15) gives the exact result, namely Eq/N = —3/8. This is not surprising 
because |3>o) = \D) at this point. At J = (Heisenberg), one has Eq/N = —0.4498, whereas the exact 
result by the Bcthe ansatz [14] is —0.4432. But one should be careful here because J = corresponds 
to one of the terminating points at which the spin-wave theory is most unreliable, as discussed below. 

From the triplet spectrum of Eq. (3.14) one sees that there are two terminating points, J = and 
1, beyond which (i.e., J < and J > 1) the spin-wave excitations are unstable. In Fig. 3, the triplet 
excitation spectrum is schematically shown for several values of the coupling constant J. It is clear that 
in the region of < J < 1, there is a nonzero gap, and this gap collapses at the terminating points. In 
particular, the triplet spectrum is flat with the gap of 1 at J = 1/2. The flatness reflects the fact that 
at J — 1/2, Hq of Eq. (3.8) contains no coupling at all between pairs of spins (dimers). More realistic 
calculations for the excitations at J — 1/2 were performed by Shastry and Sutherland and others [21]. 
They obtained the spectrum of soliton-like excitations with the minimum gap of 0.25 at k = and ir 
and the maximum gap of 1 at k = ir/2. Tonegawa and Harada's numerical calculations [19] confirmed 
the nonzero gap at J = 1/2 and in the nearby region. They predicted that the gap collapses at J w 0.3, 
while Haldane [22], who used a fermion representation, predicted this value to be about 1/6. This 
gapless point seems to signal a phase transition from the dimerized phase to a critical phase similar to 
the Heisenberg model at J = 0. In any case, the triplet spectrum of \/2sinfc from Eq. (3.14) at J = 
seems to agree well with the exact result of -| sinfc, as pointed out by Parkinson [7]. 




(3.17) 



6 



A more intriguing situation occurs for J > 1/2, where the spin- wave spectrum has a minimum at 
k = 7r/2. In particular, at J = 1, the spectrum is gapless with a cusp at k = n/2. Whether or not this 
suggests a phase change in the spatial periodicity of the system from double to four-fold, for example, is 
still unclear. The numerical calculations of the structure factor by Tonegawa and Harada [19] certainly 
showed a complicated feature for J > 1/2. There is also numerical evidence in the excitation spectra, 
which suggests that the spatial periodicity is not two- fold in the region near J = 1 [20]. Clearly, the 
spin-wave theory described here is not adequate for this task and higher-order calculations are needed. 

One can also straightforwardly calculate the long-range dimerization order parameter within the 
present spin- wave theory. The dimcr order parameter D is defined as 

D = (s,_i • b,) - (s, • b,+i) = (si(r) • s 2 (r)) - (s 2 (r) ■ si(r + 1)), (3.18) 

where the expectation is with respect to the ground state of the system. By Eqs. (2.5) and (2.6), using 
the spin- wave ground state |$o) of Eq. (3.17), it is easy to show that in the spin- wave approximation, D 
is nonzero in the region < J < 1 and gradually diminishes when J moves toward the two terminating 
points. But at the terminating points (J = 0, 1), D diverges to — oo, implying a breakdown of the spin- 
wave theory. In the following section, I provide an alternative approach to the dimerization problem by 
applying the microscopic CCM. 

IV. The Coupled-Cluster Approach 

The CCM is widely recognized nowadays as providing one of the most universally applicable, most 
powerful and most accurate of all microscopic formulations of quantum many-body theory [13]. The 
recent application of the CCM to various spin models has produced excellent numerical results [14,15]. 
It therefore seems very appropriate and timely to apply the CCM to the dimerization problem. The 
interested reader is referred to Ref. [13] for the general formalism of the CCM and to Ref. [14] for its 
particular application to the spin systems with an anticipated Ising-like long-range order. 

Generally speaking, the CCM starts with a proper model state |<j>), which is usually a simple, uncor- 
rected many-body wavefunction, and incorporates many-body correlations on top of |3>) by acting on it 
with an exponentiated correlation operator S. This operator S consists of purely so-called configuration 
creation operators with respect to the model state 1$), and is partitioned by one-body, two-body, up 
to A-body correlation operators with A the number of particles in the system. Thus, the CCM ansatz 
for the ground ket state is 

l* 9 )-c s |$). (4.1) 
The Schrodinger equation of the ground state, after a simple manipulation, can then be written as 

e - s iJe 5 |$) =E g \$), (4.2) 

where E g is the ground-state energy, and where the similarity-transformed Hamiltonian can be expressed 
as a series of nested commutators, 

e- s He s = H+[H,S\ + ± [[H, S] , S] + ■ ■ ■ , (4.3) 

which usually terminates at the fourth-order for most Hamiltonians with pair-interaction potentials [13]. 
Eqs. (4.1)-(4.3) arc the hallmarks of the CCM. 

For the case of dimerization under consideration, it is natural to choose the dimer state \D) of 
Eq. (3.2) as the model state, namely |$) = \D). The configuration creation operators with respect to 
this model state \D) are clearly given by any combinations of the three operators A[ , A r 20l and A^ . 
Their hermitian conjugates are the corresponding destruction operators. Since the antiferromagnetic 
ground state is definitely in the sector of zero s^ ot£ j , the correlation operator S is in general written in 
the form 

N/2 

S=^S n , (4.4) 

n=l 
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where N/2 is the total number of valence bonds of the spin-i chain, and the n-body correlation operators 
S n are given respectively by 



N/2 

S\ = S r A2 , 

r=l 
N/2 



O rr , Ji 10 Ji 30 2| °r,r' /i 20 /i 20 



1 



(2) 



N/2 



^(!) jr' <r" _L <?( 2 ) /tr' <r" 

-Y,r\r" /1 10 /1 30 /1 20 ~~ gj °r,r',r" yl 20^ 1 20 /1 20 



(4.5a) 
(4.56) 
(4.5c) 



etc. In Eq. (4.5b) and (4.5c) the primes on the summations imply exclusion of the terms with any pair 
of indices being equal. We notice the similarity between the spin-wave ground state |<&o) of Eq. (3.17) 
and the CCM state \^ g ) of Eq. (4.1) if S is replaced by S^- W also notice the similarity between the 
present CCM analysis and that of Ref. 15 where the spin-1 model state is given by the simple planar 
configuration with sf — for all sites I. 

The ground-state energy is obtained by taking the inner product of the Schrodinger equation (4.2) 
with the model state \D) itself, namely 

E g = (D\c- s Hc s \D); (4.6) 



the correlation coefficients {S r y,...} are determined by the coupled set of equations obtained by taking 
inner products of Eq. (4.2) with states constructed from the corresponding destruction operators, namely 

(D\A r m c- s Hc s \D) = 0, for all r, (4.7) 



for the one-body equation; 

(D\A , 01 A r 03 c- s Hc s \D) = 0, for all r,r'(^r) (4.8) 



and 

(D\A , 02 A r m c- s Hc s \D) = 0, for all r, r'& r) (4.9) 

for the two-body equations. The three-body equations and higher-order many-body equations are ob- 
tained in a similar fashion. 

One sees that the similarity-transformed Hamiltonian of Eq. (4.3) is needed in all of the above 
equations. I leave details of the derivation to the Appendix, and only point out that, as expected, the 
otherwise infinite expansion series of Eq. (4.3) indeed terminates at the fourth-order term. In fact, the 
exact energy equation (4.6) can be straightforwardly derived as 

^ = I[(l-2J)(26«+6^-a)-3], (4.10) 

where I have used the translational and reflectional symmetries, setting accordingly, 

S r = a, S^ r2 =S^ ri =b®, with i = 1,2 and r = r 2 - n, (4.11) 

It is also obvious that b^ r — bfi . 

The exact one-body equation of Eq. (4.7) can also be easily derived. It couples only to the two-body 
coefficients. Similarly, the two-body equations of Eqs. (4.8) and (4.9) couple only to the one-body and 
the three-body coefficients, and so on. From the one-body equation, it is interesting to note that the 
physical solution is given by a = 0, implying no one-body correlations for the dimerization problem. This 
is not surprising because the model state \D) is in the sector of zero total spin vector (i.e., s to tai = 0), and 
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the one-body correlation operator Si of Eq. (4.5a) will take the state out of this sector. Furthermore, if 
one assumes that the two sets of the two-body correlation coefficients are identical, namely 

b^=b?\ (4.12) 

the two-body correlation operator 5*2 then commmutes with the total spin vector s to tai- This is a 
necessary condition if one requires the CCM ground state to be in the sector of zero s to tai- (Actually, 
one requires every correlation operator S n to commute with s to tai-) The energy equation is now reduced 
to 

^ = |[(1-2J)6 1 -1]. (4.13) 

One clearly needs to employ an approximation scheme for any practical calculation. The most 
common approximation scheme in the CCM is the so-called SUBn scheme, in which one keeps up to n- 
body correlation operators and sets all the higher-order many-body correlation operators S m (m > n) to 
zero. I consider the SUB2 scheme here. I find that within the SUB2 scheme, the condition of Eq. (4.12) 
is indeed satisfied. After simplification, the two identical two-body equations are given by 

i (K3S rp + K 2 b r - 2K 1 b r+p + K X J2 K'br+p-r') = 0, r + (4.14) 

p=±l r'/O 

and 

#1 = 1-2.7, K 2 = 4(1 - 2#i6i), # 3 = #i(l + 46?)-2(l + 2J)6i. (4.15) 

A simpler approximation can be made from the full SUB2 equation (4.14), namely the so-called 
SUB2-2 scheme in which one keeps only the single coefficient, 61, setting all other b r = (|r| > 1). 
Eq. (4.14) then reduces to 

1 - 2 J + 2(3 - 2 J)bi - 9(1 - 2 J)b\ = 0, (4.16) 



with the physical solution 



h = Q(1 z 2J ) [3 ~ 2J ~ ^ 40J2 - 48J + 18] - (4 - 17) 



The full SUB2 equation (4.14) can also be solved by a Fourier transformation exactly similar to 
Eq. (3.9) of Ref. [15]. Here I only quote the final result given by the following self-consistency equation 
for bi, 

h = 3K~ 1 { 2 ~^2 1 2^ J ^Vl-^i cos2g + fc 2 cos 2 2qj , (4.18) 
where the constants k\ and k 2 arc defined by 



fci = 172 -(4#i# 2 + 8Kfh - AKfX), k 2 = U ' ^, (4.19a) 



and where X is defined by 



N/2 

X = ^ r b r b r -\.\ 

r=l 

- i /" M 



2ir 4K( cos 2q 

x (2Kx cos 2q - K 2 + K 2 y/l - k\ cos 2q + k 2 cos 2 2qj 2 . (4.196) 

After obtaining 6i as a function of J from Eq. (4.17) or (4.18), the ground-state energy is then 
given by substituting 6i into Eq. (4.13). These ground-state energies are shown in Fig. 2, together with 
the results of the spin-wave theory and of the numerical calculations [19,20] for comparison. Similar to 
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the spin-wave theory described in Sec. Ill, at J = 1/2 (the Majumdar- Ghosh point), the exact result 
is recovered for both the SUB2-2 and full SUB 2 schemes, namely, &i = and E g /N = -3/8. At the 
Heisenberg point (J = 0), the SUB2-2 and full SUB2 schemes give E g /N = -0.4268,-0.4298 respec- 
tively, slightly higher than the exact result of —0.4432. Furthermore, it is interesting to observe that in 
the full SUB2 scheme, there are also two terminating points, — —0.4443 and = 1.591, beyond 
which, namely for J < and J > Jc 2 \ there is no real solution in Eq. (4.18). The corresponding 
energy values are —0.5172 and —0.6977 respectively. The CCM SUB2 terminating points have been 
argued to correspond to the phase transition critical points in the past [14,15]. It seems reasonable to 
consider this possibility again here. From Fig. 2, one sees that the extremely simple SUB2 scheme gives 
much better results for a wide range of the coupling constant J than the spin-wave theory does, at least 
as far as the ground-state energy is concerned. 

V. Discussion 

In this paper, I have studied the dimerization problem by a microscopic approach, employing the 
proper set of composite operators of Parkinson [7]. Two approximations schemes, namely the spin- wave 
theory and a CCM analysis, have been applied to the ID frustrated spin-i model. The ground-state 
and low- lying excited energies are obtained as functions of the coupling constant. The implications of 
possible phase transitions at the naturally arising terminating points of the solutions have been discussed. 
Another approach may be provided by a variational trial wave function of the type of Eq. (4.1), similar 
to the calculation of Sachdev for the spin-i Heisenberg model [23]. 

From the present preliminary attempt to formulate a microscopic theory for the dimerization prob- 
lem, it is clear that higher-order calculations within the present analysis are needed for both the ground 
and excited states. The very successful applications [14] of the CCM to the spin systems with an 
Ising-like long-range order seem to suggest that the CCM can also provide a systematic and poten- 
tially accurate approximation scheme for the dimerization problem. Furthermore, within the formalism 
presented in this paper, it is straightforward to extend the same analysis to both higher-order dimen- 
sionality and/or spin systems with spin quantum number greater than one half. In particular, the 2D 
spin-i Heisenberg model on the square lattice with J1-J2 couplings [9,11] has been under intensive study 
for its possible dimerization. An equally interesting Hamiltonian model is provided by the ID spin-1 
Heisenberg-biquadratic systems, where it is known the ground state is dimerized at a particular coupling 
constant [4], and where trimerization is also possible in another region [10]. 
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Appendix 

In this appendix, I derive the similarity transformations within the SUB1 and SUB2 schemes of the 
CCM, described in Sec. IV. 

Notice that in the similarity transformation of Eq. (4.3), any quadratic term in the Hamiltonian of 
Eqs. (3.4)-(3.5) can be transformed as 

e-s^+V = (e- s ^.e s )(e- s A^V), (Al) 

and each similarity-transformed operator can be expanded as a series of nested-commutators, 

e- s A^e s = AT. + [A^S] + ^[[A^,S], S] + • • • . (A2) 

Since the correlation operator S consists of only the creation operators A\q, A30 and A20, the expansion 
series of Eq. (A. 2) terminates at most at the second-order by the pseudo-spin algebra of Eq. (2.4). 
Therefore the similarity-transformed Hamiltonian of Eq. (4.3) terminates at the fourth-order. 
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In the SUB1 scheme, one replaces S — > Si, where Si is given by Eq. (4.5a). From Eq. (A. 2), it is 
straightforward to derive the following SUB1 similarity transformations 



at _ Ar AT — at 

^-nO ^-riOJ "nra ^nmi 



^■oo — ^oo S r A r Q , A r n2 — A r n2 + S r A r nQ , (A3) 

^02 = ^02 + St(Aqo — A22) — S^A^q, 

where n = 1, 2, 3 and m = 1, 3, and where the definition 

ir.=e- Sl ^.e Sl , (A4) 

is used. 

In the SUB2 scheme, one replaces S — > Si +S , 2- One can firstly make the SUB1 similarity transform 
by Eq. (A. 3), and then apply the following SUB2 similarity transform for each A\-, 

^riO = ^n0' (A. 5a) 

■^nl = -^nl + 51 ^r,r'^n0^30> (A5&) 
t' 

A n2 = ^n2 - ^2 ^r.r'-^nO^O) (A5c) 
r' 

^nl = ^ri3 + *^r',T-^n0^10! (A5d) 

r' 

^00 = ^00 — (^r,r' ^10^30 + ^r',r^l0^30 — ^r,r'^20^2o)' (A5e) 
r' 

^01 = A 01 + O^r.r' ^30 (^00 — ^11 ) — ^r',r^l0^31 + ^r,r' ^20^2l) 

r' 

I 1 V^' o(!) yir' / <j(2) at 4' , "_Oc( 1 ) 4 1- at" _ c(!) 4^" \ 
"I" ^ 2^ O r,r' /i 30^°r,r" /1 20 /1 20 ZO r,r" /i 10 /1 30 °r" ,r A 10 A 30 ) 

-IT! &t S Ty' A io A lo< ~ (A5.f ) 

^02 = ^02 — O^r,r"^2o("^00 — ^22) + "Sr'.fAlO^M + AJqA^) 

t' 

— — V^' c ^ 2 - ) 4 r ' f9<?( 2 ) 4 r 4 r " — 4 r 4 r " — A r " A r \ 

2 2^ °r,r' /i 2(U z '-Vr'' /1 20 /1 20 °r,r" /1 10 /1 30 °r" ,7-^10 ^30 J 

A03 = Aj3 + E (^r',r^lo(^00 — Aj 3 ) — ^r,r'^30^13 + $r,r' ^20^23) 
1 1 V^' 47-' (c( 2 ) 4r 4 J, "_oc( 1 ) 4 r 4 r " cW 4 r " 4M 

"I" ^ 2^, °r' ,r A W\°r,r" A 20 A 20 ZO r",r /1 30 /1 10 °r,r" A 30 A W) 
r' ,r" 

where n = 1,2,3, and the primes on the summations imply exclusion of any pair of indices being equal, 
and where the operators with a bar represents the similarity transform for the S2 correlation operator, 

AI^e-^A^. (A6) 
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In deriving Eqs. (A. 5), I have used the fact that 



S$=S%=0, S™=S™. (A7) 

After the similarity transformations of Eqs. (A. 3) and (A. 5), the CCM equations of Eqs. (4.6)-(4.9) can 
be derived by using the pseudo-spin algebra of Eq. (2.4) to move in each term all the creation operators 
^10,^20,^30 to the left, all the detruction operators A i,Ao2,A 3 (also A nm , n,m ^ 0) to the right, 
and by using the fact that 

(£>K = A r mn \D) = 0, n^0, (A8) 

and 

(D\A r 00 = (D\, A r Q0 \D) = \D) 1 (A.9) 

for all r. 
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Figure Captions 



1 The perfect dimer state and dimer indexing. Each bond represents a singlet configuration as given 
by Eq. (2.1). 

2 Ground-state energy per spin as a function of the coupling constant J. Shown are results from 
the spin- wave theory (dotted), the SUB2-2 scheme (short dashed), and the full SUB2 scheme (long 
dashed). The terminating points of the SUB2 scheme and spin- wave theory are indicated. The 
numerical results from Ref. [19,20] are also included (solid). 

3 Schematic plots of the triplet excitation spectrum of Eq. (3.14) for various values of the coupling 
constant J. 
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